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I. INTRODUCTION 

There is great interest in the properties of molecules 
formed in laser-cooled atomic gases by processes such 
as photoassociation and magnetic association^^ Di- 
atomic molecules have been formed from bosonic and/or 
fermionic species of all the alkali metals from Li to Cs. 
Dimers of bosonic species have usually been found to 
decay quickly because of inelastic collisions with other 
atoms^^i but for fermionic species ( 6 Li, 40 K) it has 
been possible to stabilise the dimers by tuning to large 
positive values of the atom-atom scattering length 
This approach was used in late 2003 to create long- 
lived Bose-Einstein condensates of fermion dimers , 11 ' 12 ! 13 
and since that time there has been a large amount of 
work on their properties, and particularly on the transi- 
tion between Bose-Einstein condensation (BEC) and the 
Bardeen-Cooper-Schrieffer (BCS) regime characterised 
by long-range Cooper pairing and superfluidity. 

Photoassociation and magnetic association (tuning 
through zero-energy Feshbach resonances) both produce 
molecules that are initially in very highly excited vibra- 
tional states. However, there is great interest in bring- 
ing these molecules down to low vibrational states and 
ultimately to the absolute ground state. For example, 
Staanum et alJ^ have produced CS2 in vibrational lev- 
els v — 4 to 6 of the a 3 E + state and studied collision 
processes, while Sage et al^ have produced small num- 
bers of ultracold polar RbCs molecules in their absolute 
ground state by a 4-photon process. 
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Ultracold molecules are initially formed in the pres- 
ence of ultracold atoms and can collide with them^ For 
molecules in vibrationally excited states, there is the pos- 
sibility of vibrationally inelastic collisions, 

M 2 (v) + M — ► M 2 (V <v) + M, (1) 

where v is the vibrational quantum number. Since the 
trap depth is usually much less than 1 K, such collisions 
always release enough kinetic energy to eject both colli- 
sion partners from the trap. If the molecular density is 
high, there is also the possibility of inelastic molecule- 
molecule collisions, 

M 2 (» + M 2 (» — ► M 2 (V <v) + M 2 {v" < v). (2) 

Molecules are not destroyed in inelastic collisions, but 
they are lost from the trap and are no longer ultracold. 

The rates of inelastic atom-molecule col- 
lisions involving alkali metal dimers have 
been studied both experimentall y 14 ' 17 ' 18 ! 19 and 
theoretically^ 21 .^^^^^ Wynar et al^L formed 
87 Rb2 molecules in the second-to-last vibrational level of 
the ground excited state by stimulated Raman adiabatic 
passage (STIRAP). They estimated an upper bound of 
k\ oss = 8 x 10~ n cm 3 s" 1 due to inelastic atom-molecule 
collisions. Mukaiyama et a/. 18 measured the trap loss 
rate for 23 Na 2 molecules formed by Feshbach resonance 
tuning and obtained an atom-molecule rate coefficient 
k\ oss = 5.1 x 10~ n cm 3 s^ 1 for molecules in the highest 
vibrational state. Staanum et als^ investigated inelas- 
tic collisions of rovibrationally excited Cs 2 ( 3 S+) in 
collisions with Cs atoms in two different ranges of the 
vibrational quantum number v by monitoring trap loss 
of Cs 2 . They obtained atom- molecule rate coefficients 
close to 1.0 x 10~ 10 cm 3 s _1 for both v = 4 to 6 and 
v = 32 to 47. Zahzam et aZ. 19 carried out similar work 
for different rovibrational states of 3 S+, and also consid- 
ered molecules in the 1 S+ state and molecule- molecule 
collisions. They obtained rate coefficients of 2.6 x 10~ u 
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cm 3 s _1 and 1.0 x 10~ n cm 3 s -1 in the atom-atom and 
atom-molecule cases respectively, both with quite large 
error bounds. 

Quantum dynamics calculations on alkali metal atom- 
diatom collisions were first carried out by Soldan et al.^ 
using the potential energy surface of Higgins et al£Z- and 
subsequently extended to K + K 2 on a new potential 
surfaced In parallel work, Petrov et aZi 22 ' 23 analysed 
the stability of fermion dimers in terms of the long-range 
form of the wavefunction. In the case where the atom- 
atom scattering length a is much larger than the range 
of the atom-atom potential r e , they showed that both 
atom-molecule and molecule-molecule inelastic collision 
rates are suppressed by Fermi statistics. However, their 
derivation applies only to molecules that are in long- 
range states, with a wavefunction that depends on the 
atom-atom scattering length a, with x(r) ~ exp(— r/a). 
In the present paper, we show computationally that there 
is no systematic suppression of the atom-molecule inelas- 
tic rate for 6 Li dimers in low-lying vibrational levels, even 
when a is large and positive. We also consider mixed- 
isotope Li + Li2 collisions. This has been reported briefly 
in previous work ; 25 ' 26 but the present paper gives full de- 
tails of the calculations and further details of the results. 

The structure of the paper is as follows. Section II 
describes calculations of the potential energy surface for 
quartet Li3 , including details of both the electronic struc- 
ture methods employed and the procedures used to in- 
terpolate between and extrapolate beyond the ab initio 
points. The surface allows barrierless atom exchange 
reactions. Section III describes atom + diatom reac- 
tive scattering calculations both for homonuclear sys- 
tems [bosonic, 7 Li + 7 Li2 and fermionic, 6 Li + 6 Li2] 
and mixed-isotope systems. For collisions of 7 Li with 
either 6 Li2 or 6 Li 7 Li, exoergic atom exchange reactions 
are possible even for ground-state molecules because of 
the differences in zero-point energies. Section III also ex- 
plores the sensitivity of the cross sections to the potential 
energy surface and shows that for vibrationally excited 
states the dependence is relatively weak (only a factor of 
2 for molecules initially in v = 3). Section IV presents 
conclusions and prospects for future work. 



II. QUARTET POTENTIAL ENERGY 
SURFACES FOR LI 3 

In previous work, we investigated nonadditive forces 
in spin-polarized (quartet) alkali metal trimers 2 ^ at lin- 
ear and equilateral geometries. For Li3, we found very 
large nonadditive forces at equilateral geometries that 
reduced the interatomic distance by more than 1 A and 
increased the well depth by a factor of 4 with respect 
to the sum of Li-Li pair potentials. We subsequently 
gave a brief description of the complete potential energy 
surface 25 and identified a seam of conical intersections at 
collinear geometries between 4 E and 4 LT states. The con- 
ical intersection has been investigated further by Brue et 



al2^ The seam results in a cusp at energies close to the 
atom-diatom threshold and influences Li + Li2 collision 
dynamics at ultralow energiesi 25 ' 26 In the present paper, 
we give a full description of the ab initio calculations and 
fitting procedure used to obtain a global representation 
of the surface and extend the collision calculations. 



A. Electronic states overview 

In the present work, we are interested mostly in the 
lowest quartet electronic state of Li3, which is a 4 A' 
state in C s symmetry. This state correlates with three 
ground-state ( 2 5 l ) atoms at large interatomic separations. 
As mentioned above, the quartet ground state has cross- 
ings with other electronic states that correlate with the 
S + S + P separated-atom limit. These crossings (and 
avoided crossings) occur at energies that are potentially 
relevant to scattering processes involving three ground- 
state Li atoms. To gain a broader, qualitative picture on 
how the different potential energy surfaces intertwine, 
we have carried out state-averaged multiconfiguration 
self-consistent field (MCSCF) calculations of the poten- 
tial energy surfaces of quartet states for certain high- 
symmetry arrangements. These calculations used the 
complete active space (CASSCF) method of Werner and 
Knowles^ii 3 ^^^ as implemented in the MOLPRO 
package^ with an aug-cc-pVTZ basis set generated from 
the cc-pVTZ set 3 ^ as described below. The active space 
included all 24 molecular orbitals constructed principally 
from 2s, 2p, 3s and 3p atomic orbitals, with the three Is 
orbitals doubly occupied in all configurations; this active 
space is designated [3,24] below. 

Figure [T] shows the potential curves obtained for quar- 
tet states at Dooh geometries. The ground state is of E u 
symmetry; at a Li-Li bond length r near 3 A it intersects 
a Tig state, which in turn intersects a E ff state at higher 
energies. If we change the interbond angle a, keeping 
the two adjacent bond lengths unchanged (moving from 
Dooh to Civ arrangements), the ground state becomes 
of £?2 symmetry, while the U g state splits into states of 
A2 and B2 symmetries. The two B 2 states avoided-cross 
around a — 70°. At D^h arrangements (a — 60°), shown 
in Figure [21 the ground state is of A' 2 symmetry, while 
the second B 2 state forms a doubly degenerate E' state 
with the A\ state that originated from S ff at -Doo/i- The 
A2 state in C% v symmetry (LI g = A 2 + B 2 ) forms an E" 
state with the B\ state that was a part of H u = A\ + B-y 
at -Doo/i- 

For comparison, Figure [3] shows the triplet states of 
Li 2 calculated using [2,8]-CASSCF, followed by an inter- 
nally contracted multireference configuration interaction 
(MRCI) 34 i 37 ' 38 i 39 including single and double excitations 
from the CASSCF wavefunction. 

Figure [4] shows a correlation diagram that connects 
the Dooh, D^h and atom-diatom {Civ) limits for states 
of quartet Li3 that correlate with the S + S + P atomic 
limit for fixed bond length of 6 A. Note that the molecular 
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FIG. 1: CASSCF quartet potentials of L13 at geometries 
for states that correlate with the atomic S+S+S and S+S+P 
asymptotic limits. The interatomic distances are r, r, and 2r. 



FIG. 3: Triplet potential energy curves of Li2 from the atomic 
S + S and S + P dissociation limits. 
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FIG. 2: CASSCF quartet potentials of Li3 at D^u geometries 
for states that correlate with the atomic S+S+S and S+S+P 
asymptotic limits. The interatomic distances are r. 



state arising from three S'-state atoms is not shown. 

In the present work we focus on calculating the lowest 
quartet state, 4 A' , to high precision. Colavecchia et a/4°- 
and Brue et al£^ have carried out full configuration inter- 
action (FCI) calculations with just one electron on each 
Li atom correlated. Our calculations differ from theirs in 
correlating all electrons and in using considerably larger 
basis sets, but use approximate correlation treatments. 



B. Basis set convergence 

A variety of basis sets are available for lithium. Feller 36 
constructed correlation-consistent polarized valence ba- 
sis sets (cc-pVXZ for X=D, T, Q, and 5)^1 and Iron et 
al£^ devised correlation-consistent polarized core- valence 
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FIG. 4: Correlation diagram of quartet potentials of Li3 that 
correlate with the atomic S + S + P dissociation limit. Note 
the the energy is shown relative to this limit and the state 
correlating with S + S + S is not shown. The first panel 
connects the -Dooh states with the atom-diatom limit, with 
one interatomic distance fixed at 6 A and C xv symmetry 
preserved. The second panel connects the atom-diatom limit 
with Dzh states, with one interatomic distance fixed at 6 A 
and C<2 V symmetry preserved. The third panel connects D^h 
and Daoh terms, with two interatomic distances fixed at 6 A 
and the angle between them, a, varied. 



basis sets (cc-pCVXZ) by adding "tight" functions to 
Dunning-type cc-pVXZ bases in order to give a better 
account of the core-core and core-valence correlation. 
There are no standard augmented basis sets available 
for Li, — so we generated augmented variants of several 
basis sets using the even-tempered scheme implemented 
in MOLPRO, which adds an additional diffuse function 
of each angular symmetry (s, p, d, f, ...), with an even- 
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tempered exponent based on the ratio of the two smallest 
exponents in the original set (or a ratio of 2.5 if only one 
function of a type is present in the original basis). 

The cc-pVXZ and cc-pCVXZ basis sets are usually 
used with the s and p functions contracted. However, 
in previous work on alkali metal trimers^ we found it 
advantageous to use the basis sets in uncontracted form. 
In the present work we consider the cc-pV5Z basis set in 
three different forms: fully contracted (sp-contracted), 
fully uncontracted, and with just the p space uncon- 
tracted (s-contracted). However, we found it too expen- 
sive to uncontract the cc-pCVXZ and augmented basis 
sets and these were used in fully contracted form. 

We chose to use RCCSD(T) calculations (restricted 
open-shell coupled cluster with single, double and nonit- 
erative triple excitations). All calculations were carried 
out using MOLPRO^ 

With the computer resources available to us, we found 
that RCCSD(T) calculations of triatomic energies with 
the largest basis sets available, cc-pCV5Z and aug-cc- 
pCV5Z, were feasible only for relatively small numbers 
of points, not for a grid of several hundred points. These 
basis sets were therefore used only for atomic and di- 
atomic calculations, except that the cc-pCV5Z basis set 
was used to provide benchmark triatomic calculations at 
a restricted number of points. 

We began by testing convergence of two atomic prop- 
erties, the static polarizability and the S—P excitation 
energy, which are important for accurate calculations of 
long-range forces. The results for a variety of basis are 
shown in Table |TJ In general terms all the basis sets 
perform acceptably for both quantities, though the s- 
contracted cc-pV5Z basis set is fortuitously accurate for 
the S—P excitation energy. 

We next considered the performance of different basis 
sets on potential curves for the 3 E+ state of L12. Elec- 
tronic energies were calculated on a mesh of 64 inter- 
atomic distances. The results were corrected for basis- 
set superposition error (BSSE) using the counterpoise 
correction^ The points were then interpolated using the 
reciprocal power reproducing kernel Hilbert space (RP- 
RKHS) interpolation metho d 47 ! 48 with dispersion coef- 
ficients C§, Cg and C10 fixed to accurate values from 
the calculations of Yan et al^ At internuclear distances 
r > 16 A the ab initio energies were replaced by those 
from the three-term dispersion formula. Other interpola- 
tion parameters, in the notation of Ref. l48l . were m = 2, 
n = 3, r a = 15 A, r 65 = 21.5 A, r 66 = 22.5 A, r 67 = 23.5 

A. 

The resulting diatomic well depth D e and bond length 
r e are included in Table HI together with the binding en- 
ergy of the last bound state E\q and scattering length a 
calculated for lA?. The results are compared with prop- 
erties obtained from the RKR curves of Linton et al*& 
(the values of a and E\q are from photoassociation spec- 
troscopy in combination with RKR results 51 ). It may 
be seen that all the basis sets underestimate the RKR 
well depth by between 2 and 5 cm -1 , except for the cc- 



pCVQZ basis set which underestimates it by 8 cm -1 and 
the sp-contracted cc-pV5Z basis set which slightly over- 
estimates it. On this basis we excluded the cc-pCVQZ 
basis set from further consideration. 

Another way to test the quality of a basis set is to con- 
sider the magnitude of BSSE. BSSE arises when basis 
functions on two centers overlap and compensate for in- 
adequacies in the one-centre basis. Table U shows dimcr 
BSSE values at r — 4.2 A, which is near the potential 
minimum. The cc-pCVXZ basis sets give quite small 
BSSE values since they are by construction geared to 
represent the core-core correlation well. However, as dis- 
cussed above, the cc-pCVQZ basis set is not accurate 
enough for our purposes and the cc-pCV5Z basis set is 
too expensive. We therefore focussed on the cc-pV5Z ba- 
sis sets. As may be seen in Table U the cc-pV5Z basis 
set in sp-contracted form gives a very large BSSE, 41.45 
cm -1 at r = 4.2 A, which we considered unacceptably 
high. However, uncontracting the p space reduces this to 
2.553 cm -1 . Uncontracting the s space as well reduces 
this by a further factor of 3, but at a very considerable 
computational cost. 

Finally, we carried out convergence tests on triatomic 
energies. In this case there are no experimental results to 
compare with, so we used the position and energy of the 
equilateral minimum, calculated with the cc-pCV5Z basis 
set, as our benchmark. As may be seen in Table[IJ there is 
very little difference between the s-contracted and uncon- 
tracted forms of the cc-pV5Z basis set, which underesti- 
mate the cc-pCV5Z well depth by 11 and 13 cm" 1 respec- 
tively. However, the sp-contracted form overestimates 
the well depth by 40 cm" 1 . For triatomic RCCSD(T) cal- 
culations the s-contracted basis set costs only 10% more 
than the sp-contracted set, which the fully uncontracted 
set costs twice as much. On this basis we decided to pro- 
ceed with the s-contracted cc-pV5Z basis set to calculate 
the full potential surface. 

One further possibility we considered was the use of 
frozen-core calculations, which are considerably cheaper 
than calculations that correlate all electrons. However, 
we found that such calculations systematically overesti- 
mate dimer well depths by 2 to 3% and equilibrium dis- 
tances by about 2% compared to cc-pCVXZ calculations 
at each level (T, Q, Z). For example, frozen-core cal- 
culations with the sp-contracted cc-pV5Z basis set give 
D e = 338.471 cm" 1 and r e = 4.199 A for Li 2 . Frozen- 
core calculations on the trimer with this basis set un- 
derestimate the trimer well depth, giving D e = 3937.5 
cm" 1 and r e = 3.137 A. We therefore did not pursue the 
frozen-core approximation. 

As described above, the Li3 potential shows a seam 
of conical intersections at linear geometries. The 
RCCSD(T) approach is inherently a single-reference 
method, so its appropriateness in the vicinity of a conical 
intersection is open to question. The usual approach is to 
consider the Tl diagnostic, which can be used to identify 
where multireference effects become large and may com- 
promise the results. Lee and Taylor—^ have given a rule of 
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Atomic 




Diatomic 










Triatomic 




basis set 


a(og) 


S—P energy 


re (A) 


D e 


a (ao) 




BSSE 


re (A) 


D e 






(cm" 1 ) 




(cm" 1 ) 




(cm 1 ) (cm 1 ) 




(cm" 1 ) 


sp-contracted cc-pV5Z 


165.518 


14834.01 


4.169 


334.042 






41.45 


3.102 


4021.5 


s-contracted cc-pV5Z 


164.396 


14903.62 


4.177 


328.922 


12.54 


0.223 


2.553 


3.102 


3969.8 


uncontracted cc-pV5Z 


164.336 


14906.69 


4.177 


328.732 






0.8015 


3.102 


3967.4 


aug-cc-pCV5Z 


164.189 




4.176 


330.548 


-8.95 


0.346 


0.234 






cc-pCV5Z 


164.152 


14906.48 


4.175 


328.952 






0.139 


3.098 


3980.6 


cc-pCVQZ 


164.140 


14911.28 


4.179 


325.499 






0.501 






Best available theory 


164.111" 


















Experiment 


164.0±3.4 b 


14903.89 c 


4.173 


333.69 


-27.3 ± 0.8 


0.416 









TABLE I: Convergence tests for Li basis sets. For atoms: static polarizability a and S—P excitation energy. For triplet dimers: 
dissociation energy D e , position of the minimum r e , BSSE evaluated at r = 4.2 A scattering length a and energy of the highest 
vibrational level Eio for the 7 Li2 molecule. For quartet trimers, dissociation energy D e and position of the minimum r e for Li3 
(Dih)- All calculations except "best available theory" are from the present work and used RCCSD(T) calculations. a : ref. |43|; 
b : ref. H c : ref. EE 



thumb that single-reference methods become unreliable 
when Tl > 0.02. In the present calculations, the Tl diag- 
nostic was relatively small at near-linear geometries but 
had a maximum value of around 0.11 at equilateral ge- 
ometries. We therefore investigated the reliability of the 
RCCSD(T) approach by recalculated the Li 3 potential 
energy surface with a smaller basis set (aug-cc-pVTZ), 
using [3,12]-CASSCF followed by a multireference config- 
uration interaction (MRCI) calculation including David- 
son's correction, and compared it with the corresponding 
RCCSD(T) energies. The Tl diagnostic was again large 
in the well region for bond angles of ~ 60°. Nevertheless, 
the two surfaces agreed within 2% in well depth with no 
visible qualitative differences. We therefore believe that 
the RCCSD(T) energies obtained with the cc-pV5Z basis 
set are trustworthy. 

C. The 1 4 A' surface 

1. Representation of the surface 

If the potential is to be capable of representing all the 
properties of experimental interest (including atom-atom 
scattering lengths, dimer and trimer bound states, atom- 
diatom collisions and 3-body recombination), then it is 
very important that it should represent dissociation into 
all possible sets of products (atom + diatom and 3 sep- 
arated atoms) with the correct long-range behaviour. In 
our previous work on K + K2, 24 this was achieved by 
separating the potential into pairwise-additive and non- 
additive parts and interpolating them separately, 

3 

V(r) = J2Vdim e r(r ij )+V 3 (r), (3) 

i<j 

where r represents (ri 2 , r-iz , rax ) and the functional form 
for the nonadditive part V^r) is chosen to ensure the 



correct behaviour in the atom-diatom limits However, 
for Li + Li 2 the non-additive term is so large that the 
potential minimum for the trimer occurs at a distance 
that is high on the repulsive wall for the dimer. Because 
of this, Eq. ^ would represent the interaction poten- 
tial in this region as a difference between two very large 
quantities. This is unsatisfactory. Nevertheless, at long 
range a decomposition according to Eq. (J3J) is essential. 
Under these circumstances, it is best to fit the ab initio 
points directly to obtain a short-range function Vsit(r) 
without imposing the correct long-range behaviour, and 
to join this onto a potential with the proper long-range 
behaviour VLR,(r) using a switching function S(r), 

V(r) = S(r)V SR (r) + [1 - %)]V LR (r). (4) 

The switching function is 1 at short range but switches 
smoothly to zero at long range as described below. 

2. Choice of grid 

For quantum dynamics calculations, it is very impor- 
tant to represent the potential energy function smoothly 
and without oscillations between ab initio points. There 
are several coordinate systems that can be used for tri- 
atomic systems, including hyperspherical coordinates, 
Jacobi coordinates, valence (bond length / bond angle) 
coordinates and pure bond- length coordinates. These are 
by no means equivalent for interpolation purposes. As 
will be seen below, the dynamical calculations are car- 
ried out in hyperspherical coordinates. However, grids 
of points in hyperspherical coordinates tend to include 
points in which 2 atoms lie very close together, which 
hinders interpolation because polynomials with very high 
localised maxima tend to have oscillations in other re- 
gions. Jacobi coordinates suffer from the same problem, 
and neither Jacobi coordinates nor valence coordinates 
allow the full 3-body exchange symmetry to be intro- 
duced in a natural way. Fortunately, there is no need to 



represent the potential energy surface in the same coordi- 
nate system as is used in the dynamical calculations. We 
therefore chose to carry out electronic structure calcula- 
tions on a grid of points in pure bond-length coordinates 
(^12,^23,^31)7 which make exchange symmetry easy to 
handle but avoid points with atoms too close together. 

The grid was constructed using the internuclear dis- 
tances r = 2.0, 2.4, 2.8, 3.2, 3.6, 4.0, 4.4, 4.8, 5.2, 6.0, 
6.8, 8.4 and 10.0 A. At nonlinear geometries, the grid in- 
cluded all 315 unique combinations of (ri2, ^23 and f3i) 
taken from this set that satisfy the triangle inequality, 
r 3i < r i2 + r 23- We also added additional points for 
which r\i and r23 are in this set and 7-3 1 takes values 
larger than 10.0 A in steps of 1.6 A. This produced an- 
other 56 points, such as, for example, (6.0,6.0,11.6) A. 
A grid of 120 linear configurations (rai = ri2 + 7-23) was 
formed by taking all possible combinations of r\% and 7-23 
from the above set with the additional distance 5.6 A. We 
calculated electronic energies at a total of 491 geometries. 

Interaction energies for the 4 A' state were calculated 
with respect to the 3-atom S + S + S dissociation limit 
using RCCSD(T) calculations with the s-contracted cc- 
pV5Z basis set described above. All interaction en- 
ergies were corrected for BSSE using the counterpoise 
correction^ 



3. Fitting and interpolation 

Even with almost 500 points, the grid is fairly sparse 
and the quality of interpolation is important. However, 
many interpolation procedures in N dimensions work 
only for a "product" grid, in which points are available at 
all combinations of a set of coordinate values. This can- 
not be provided in pure bond length coordinates because 
of the need to satisfy the triangle condition. We therefore 
chose to use the interpolant-moving least squares (IMLS) 
method. 54 In this approach, approximate gradients and 
Hessians are calculated at each grid point and stored. 
The fitting function is represented as a weighted sum of 
Taylor expansions of second order about the grid points. 

The IMLS procedure has two parameters, the power p 
that determines the shape and range of the weight func- 
tion and a smoothing parameter e that removes the flat 
spots at data points. Figure 5 shows IMLS fits to the 
potential for D^h geometries in the region of the conical 
intersection, with p = 6 and e values of 0.03 and 0.05 A. 
The grid points and the surface of Colavecchia et al^- are 
also shown for comparison. The value of e used here is 
large in comparison with previous application s 54 ! 55 ' 56 ! 57 
and produces a root-mean-square error of 9.72 cm -1 . 
The relative error was larger than 1% at 67 geometries, 
21 of which lie below the 3-body dissociation limit and 
are enclosed by 0.5 contour of the switching function (see 
Eq. ||7J) below). There are only 5 geometries at which the 
relative error is larger than 10%, and these lie either near 
the conical intersection or near the zero of the potential 
(where the relative error does not reflect the accuracy). 
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FIG. 5: Comparison of different fitted potentials with the ab 
initio energies for the quartet ground state of Li3 at D^h 
geometries. 



4- The long-range potential 

The long-range potential, VLrtXr), was constructed as 
a sum of pairwise potentials and a nonadditive part as in 
Eq. ([3]). The need to include the nonadditive term has 
been described in detail in ref. l53l. In particular, the sum 
of pair potentials gives an isotropic atom-diatom Ce coef- 
ficient, whereas in reality there are important anisotropic 
terms. For the pairwise potentials Vdimcr(?") we used the 
Li2 RKR potential of Linton et al., 50 which extrapolates 
to a three-term dispersion expression with coefficients Cq , 
Cg and C10 as given by Yan et al^ The nonadditive po- 
tential was constructed as described in reb^ 3 - from the 
Axilrod-Teller-Muto triple-dipole termj 5 ^ 5 ^ the third- 
order dipole-dipole-quadrupole dispersion energ y 60 ! 61 
the fourth-order dipole dispersion energ y 62 ! 63 and an 
atom-diatom term that is short-range in one of the three 
distances,— 

^3,rc P (r)= - [A + BP 2 {cos6)]cM-Cx)r^r^ 
x [D 6 (r 23 )D 6 (r 31 )]^ 2 
+ cyclic permutations, (5) 

where x = (r%2 —ro)/ra and 9 is a Jacobi angle approxi- 
mated in terms of internal angles (fi as 

P 2 (cos6>) ~ — -(1 + 3cos</?i cos(y9 2 cos^). (6) 

All the non-additive terms were damped with a combi- 
nation of Tang-Toennies damping functions Dq and Dg 64 
as described in ref. l53l. Three-body dispersion coefficients 
Z\\i and Z\\2 were taken from refs. |43| and IrJEI. respec- 
tively, Zxiix was estimated as described in ref. [H, and 
the parameters A, B and C were determined as described 
in ref. 53 by fitting to the isotropic and anisotropic disper- 
sion coefficients for Li + Li 2 ( 3 S^) obtained as a function 
of r by Rerat et alf& 
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2 3 4 5 6 7 8 

^(A) 

FIG. 6: The fitted quartet ground-state potential of lithium 
trimer for a bond angle of 60°. Contours are labeled in cm - . 



The long-range form Vlr(t) is designed to be valid 
when any of the atom-atom distances is large. The short- 
range (IMLS) and long-range potentials were joined using 
the switching function 

S(r) = i tanh[l - s 1 (r 12 + r 2 i + r 31 - s 2 )], (7) 

with si set to 0.7 A -1 and s 2 to 20 A. 

Contour plots of the final potential for three different 
interbond angles a are shown in Figures [6] to [8] 

III. LI + LI 2 SCATTERING CALCULATIONS 
A. Methodology 

Both linear and equilateral configurations of the Li3 
collision complex (and all geometries in between) lie lower 
in energy than Li + Li2. The potential energy sur- 
face thus allows barrierless atom exchange (rearrange- 
ment) collisions even at limitingly low collision ener- 
gies. Calculations using an inelastic formalism (includ- 
ing only one arrangement channel) are therefore inade- 
quate. Accordingly, we have carried out quantum reac- 
tive scattering calculations to obtain collision cross sec- 
tions. The methodology has been described in our pre- 
vious work on ultracold collisions in Na + Na 2 20 i 21 and 
K + and also in studies of insertion reactions such 
as N( 2 L>) + H 2 -> NH + H, 67 0( 1 L>) + H 2 -» OH + 
and S( 2 D) + H 2 -> SH + E££ at thermal energies. 

The time-independent Schrodinger equation is solved 
using a coupled-channel method in body-fixed demo- 
cratic hyperspherical coordinates using a diabatic-by- 




2 3 4 5 6 7 8 

^(A) 

FIG. 7: The fitted quartet ground-state potential of lithium 
trimer for a bond angle of 120° . Contours are labeled in cm" 1 . 




2 3 4 5 6 7 8 

MA) 

FIG. 8: The fitted quartet ground-state potential of lithium 
trimer for a bond angle of 150° . Contours are labeled in cm" . 



sector method. In each sector, the wavefunction is ex- 
panded in the eigenstates of a reference hamiltonian 
at a fixed hyperradius. The reference hamiltonian in- 
cludes potential energy terms and kinetic energy terms 
arising from deformation and rotation around the axis 
of least moment of inertia. Its eigenfunctions are ob- 
tained by variational expansion in a set of pseudohy- 
perspherical harmonics. The resulting coupled equa- 
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tions are integrated using the log-derivative method of 
Manolopoulos^ At the point where interchannel cou- 
plings are small enough to be neglected, the coupled- 
channel solutions are matched onto atom-diatom func- 
tions propagated inwards in Jacobi coordinates to take 
account of the isotropic long-range and centrifugal po- 
tentials. 

The diatom wavefunctions for Hund's case (b), which is 
appropriate for Li2 ( 3 S+), arc labeled with a vibrational 
quantum number v and a mechanical rotational quantum 
number n; n couples with the diatomic electron spin s = 
1 to give a resultant j. However, for Li2 the splitting 
between states of the same n but different j are very 
small. We therefore do not introduce the spin explicitly 
and n and j are equivalent in our calculations. 

The matching procedure yields the reactance matrix 
K, the scattering matrix S and the transition matrix 
T = I — S, where / is a unit matrix. The differential 
cross sections, averaged over initial m n and summed over 
final m' n states, are given by 



The cross sections can be written exactly in terms of the 
energy-dependent scattering length, 



(da 

V / rvn—^T'v'r 

l 

4(2n + l)k 2 v 



JKK' 



where r labels the arrangement, v and n indicate the vi- 
brational and rotational quantum numbers of the diatom, 
K labels the z component of J and n in the body-fixed 
frame, and k Tvn is the corresponding wave vector. Q is 
an element of solid angle at a scattering angle i? with 
respect to the initial approach direction in the centre- 
of-mass frame and the d functions are reduced Wigner 
rotation matrices. 

Integral state-to-state cross sections are obtained by 
integrating Eq. [8j 



(2n + l)k 2 



vn JKK' 



(2J+1) It, 



r'v'n'K'rvriK 



(9) 

At ultralow energies, only Tu matrix elements with atom- 
diatom end-over-end angular momentum / = contribute 
significantly to the cross sections. The elastic and total 
inelastic cross sections in terms of Tu are then 



Ocla 



k 2 



\Tu\ 2 , 



'incl 



^Li-li- T « 



(10) 

(ii) 



where the rvn subscript has been dropped to simplify 
notation. At low energy, Su and Tu are conveniently 
parameterized in terms of an energy-dependent complex 
scattering length a(k) — a(k) — i(3(k), where 



f,\ 1 ( I- Si, 

a{k) = -k\TTs-, 



1 



i/c \ 2 Tj'i 



(12) 



47r|ar 



^clas 



incl 



l + k 2 \a\ 2 + 2kf3' 
4tt/3 

fc(l + k 2 \a\ 2 + 2k/3)' 



(13) 
(14) 



The scattering length becomes constant at limitingly low 
energy. If ka <C 1, which is true at limitingly low energy 
except near a zero-energy resonance, Eq. rT2] reduces to 
the expression commonly used for the zero-energy scat- 
tering length j2S 



1 T u 
a = — lim — . 

2i fc-^o k 



(15) 



Elastic and inelastic cross sections at ultralow energies 
are governed by the Wigner threshold laws,— 



(8) 



°clas ~ E 21 , 

jpl-X/2 
CTinel ~ , 



(16) 



where E is the collision energy and I is the end-over-end 
angular momentum. In the presence of a long-range po- 



tential —C. 



elastic cross sections have an additional 



term ~ E s 3 that dominates for higher I at ultralow en- 
ergies. 



B. Computational details 

We have performed coupled-channel atom-diatom col- 
lision calculations on all possible mixtures of the two iso- 
topes of lithium, 7 Li and 6 Li. The 7 Li nucleus is a fermion 
with nuclear spin 3/2, so that the atom, including its elec- 
trons, behaves as a composite boson in cold dilute gases. 
The 6 Li nucleus is a boson with nuclear spin 1, while the 
atom is a composite fermion. The present calculations 
are restricted to collisions in which all three of the atoms 
are in spin-stretched states (with F = I + S = F max and 
| | = F). For these states the nuclear spin wavefunc- 
tion is symmetric with respect to any exchange of nu- 
clei. However, the electronic wavefunction for the quar- 
tet state, which depends parametrically on nuclear co- 
ordinates, is antisymmetric with respect to exchange of 
two identical nuclei. This means that the spatial nu- 
clear wavefunction must be symmetric under exchange 
of any two fermionic nuclei (bosonic alkali metal atoms) 
and antisymmetric under exchange of any two bosonic 
nuclei (fermionic alkali metal atoms) . The symmetry un- 
der exchange of identical nuclei is easy to impose in our 
program by selecting pseudohyperspherical harmonics of 
appropriate symmetry in the basis set. 

Propagation of the coupled-channel solutions is per- 
formed separately for each partial wave J n , labeled by 
the spin-free total angular momentum J, and parity 
eigenvalue n. The parity is (— l) n+l , and only I values 
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that satisfy the triangle inequality, | J — n\ < I < J + n, 
are permitted. 

The basis set for 7 Li 3 included all eigenstates of the 
reference hamiltonian (see above) that match onto 7 Li2 
states below the v = 7 manifold. Only even values of n 
are allowed for 7 Li2 in spin-stretched states. The result- 
ing number of coupled equations varied between 97, for 
total angular momentum J — 0, and 827, for J = 10. For 
the 6 Li 3 system, including only odd-n states of 6 Li 2 be- 
low the v = 7 manifold gave a basis set that varied in size 
between 85 for J = and 782 for J = 11". For 6 Li 7 Li 2 
and Li 6 Li 2 , where the symmetry is reduced, the basis 
set even for J = consisted of 272 and 263 functions, re- 
spectively. The calculation of eigenstates of the reference 
hamiltonian involved diagonalizations of matrices of sizes 
between 1240 and 2136 for 7 Li 3 , between 1220 and 2162 
for 6 Li 3 , and between 3660 and 6488 for the isotopically 
mixed systems. 

The coupled equations were integrated across a range 
of hyperradius p between 5 eto and 45 ao with a log- 
derivative sector width of 0.1 ao. Within each sector, 
the wavefunction is propagated using the log-derivative 
method of Manolopoulos, 72 with a step size keyed to 
the local de Broglie wavelength. The propagation in 
the outer region in Jacobi coordinates extended as far 
as R = 10000 ao for collision energies of E = 1 nK. 



C. Ultracold collisions 

In the ultracold regime, the dominant contribution to 
cross sections comes from the partial wave that includes 
the I = channel (s-wave scattering). For a particular 
initial rovibrational state (v, n), the cross section is thus 
calculated from the partial wave labeled by the total an- 
gular momentum J — n and parity II = (—1)™. Figures 
l9l and fTTJl show the J = cross sections for 7 Li + 7 Li 2 for 
initial states (v,n) — (0,0) and (1,0). The insets show 
the real and imaginary part of the complex scattering 
length, calculated using Eq. (Tl2|) . 

The energy dependence is particularly simple for col- 
lisions involving a ground-state 7 Li 2 molecule, (i>, n) = 
(0, 0). Only elastic scattering is allowed below 2.3 K (the 
energy of the (0, 2) threshold). The elastic cross sections 
are essentially constant below ~ 10~ 5 K, and oscillate at 
higher energies. The s-wave cross sections exhibit sud- 
den drops to zero at 37 mK and 300 mK and a feature 
around 750 mK. The one-channel scattering is character- 
ized by a phase shift 8 that can be extracted from the 
1 x 1 K matrix element, Kn = tan The energy de- 
pendence of the phase shift is shown in Figure [Til an d 
the zeroes in the cross section can be associated with 8 
passing through a multiple of n. The profile of the cross 
section follows ~ sin 2 8. However, it should be noted that 
partial waves with I > make significant contributions 
at collision energies above about 1 mK, so that the zeroes 
in the s-wave cross section will be washed out in the total 
elastic cross section. 




FIG. 9: Energy dependence of elastic cross sections for 7 Li + 
7 Li2(vi = 0, m = 0). The energy-dependent scattering length 
is shown in the inset. 
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FIG. 10: Energy dependence of elastic and inelastic cross 
sections for 7 Li + 7 Li2(i>i = l,Hi = 0). The complex energy- 
dependent scattering length is shown in the inset. 



The feature at ~ 750 mK is different. It is associ- 
ated with a sharp rise in 8 through n, superimposed on 
the falling background. This is characteristic of a scat- 
tering resonance, and in this case is due to a Feshbach 
resonance involving a quasibound state associated with 
a rotationally excited threshold. 

For the vibrationally excited initial state (v, n) = (1, 0) 
the above arguments generalize. The S-matrix is now an 
9x9 matrix with 8 inelastic (or reactive) channels cor- 
responding to v = 0, n = 0, 2, . . . 14. The elastic cross 
sections, Figure[T2l still show dips, but these do not reach 
zero because the elastic S-matrix element cannot be 1 in 
the presence of inelastic scattering. The quantity that 
now shows a rise of n across a resonance is the eigen- 
phase sum, which is the sum of phases obtained from 
eigenvalues of the K (or S) matrix. The eigenphase sum 
and the individual eigenphases are shown in Figure 1131 



O i 




o r i i i i i i i i i i i i i i i i i i i i i i , , 1 
0.2 0.4 0.6 0.8 1 

E(K) 

FIG. 11: Energy dependence of the phase shift 8 for 7 Li + 
7 Li 2 (vi = 0,m = 0). 
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FIG. 12: Energy dependence of elastic and sum of inelastic 
probabilities, ak 2 /ir, for 7 Li + 7 Li2(vi = l,m = 0). 



It may be seen that there are two broad resonances cen- 
tered at about 0.5 K and 0.8 K, but that the phase change 
associated with them is spread between all the channels. 
The individual eigenphases show a complicated pattern 
with many avoided crossings. 

The inelastic cross sections show a slope of —1/2 on 
a log-log scale at ultralow temperatures, as predicted by 
the Wigner threshold law, Eq. (fT7J]) . Above about 10~ 2 
K, the slope changes to —1. The total inelastic prob- 
ability, shown in Figure 1121 increases with collision en- 
ergy according to the Wigner laws and saturates close 
to unity outside the Wigner region^ At this point the 
energy dependence of the cross sections is governed by 
the kinematic prefactor l/k 2 in Eq. (fTTj) . A high inelas- 
tic probability is a feature of reactions proceeding over a 

deep weu.aa^iSLsa^za 

The final state rotational distributions at ^ 1 nK are 
shown in Figure 1141 The oscillatory structure is simi- 




i i i i i i i i i i i i i i i i i_ 

0.2 0.4 0.6 0.8 

E (K) 



FIG. 13: Eigenphase sum (upper panel) and individual eigen- 
phases (lower panel) for 7 Li + 7 Li2(vi = l,m = 0). 



lar to that found earlier in ultracold collisions of Na + 
Na2^ and similar structures have been observed in vibra- 
tional predissociation of Van der Waals complexes. The 
oscillation is thought to be due to a rotational rainbow 
effect. 76 In classical terms, the angular momentum im- 
parted to the hi2 molecule is zero if the kinetic energy 
is released at linear or T-shaped geometries, but large at 
around 9 = 45°. In this model, the oscillations arise from 
interference between classical trajectories on either side 
of the maximum. 

The relative rotational distributions do not change sig- 
nificantly across the Wigner regime, although individual 
cross sections vary over five orders of magnitude between 
1 nK and 100 mK because of the dependence of the 
Wigner threshold laws. 

We have studied collisions involving a range of different 
initial molecular states. Results at a collision energy of 
~ 1 nK are presented in Table HT1 There is no systematic 
dependence of the cross sections on the initial molecular 
state. Instead, the cross sections show essentially random 
scatter about a mean value. The origin of this behavior 
is discussed below. At 1 nK, the inelastic cross sections 
are typically three orders of magnitude higher than the 
elastic cross sections. The ratio decreases with the colli- 
sion energy according to the Wigner laws, Eq. (fTS)) . until 
it reaches ~ 1 in the millikelvin range. Sympathetic and 
evaporative cooling of atoms and molecules depend on 
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FIG. 14: Final rotational distributions for 7 Li + 7 Li2(vi 
l,m = 0) at 0.928 nK. 



a low value (< 10 -2 ) of this ratio, so it is unlikely that 
such cooling mechanisms will be possible for alkali metal 
dimers in low-lying vibrationally excited states. 

It has previously been found that collisions involving 
rotationally excited diatoms can cause unusually efficient 
and specific energy transfer when the collision time is 
longer than the rotational period and there is a match 
between the rotational and vibrational frequencies i 77 > 78 
This has been termed quasiresonant energy transfer 
to distinguish it from energy-resonant behaviour aris- 
ing from energy matching between the initial and final 
states. Quasiresonant energy transfer has been explained 
in terms of adiabatic invariance theory 7 ^ with the as- 
sumption that the coupling between the degrees of free- 
dom is not strong. Such behaviour has also been stud- 
ied at ultralow temperatures i 80 i 81 ' 82 i 83 Inelastic transi- 
tions in some highly excited rotational levels of hydro- 
gen molecules were found to be dramatically suppressed 
when quasiresonant transitions were energetically forbid- 
den. However, the highest initial rotational states consid- 
ered for 7 Li2 have j = 10, for which the rotational energy 
is only about half the vibrational spacing. We thus do 
not expect to observe quasiresonant enhancements of the 
type described in refs. [7?] and In addition, we ob- 
serve no near-resonant enhancement of the inelastic rate 
for initial state v = 2, j = 0, for which a transition to 
the v = 0, j = 20 level is exothermic by only 3.33 K. As 
may be seen in Table HH the inelastic cross section from 
this level is actually the lowest of all those calculated. 

The atom-molecule inelastic loss rates observed in ex- 
periments on ultracold boson dimers ( 133 Cs2) 
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which is consistent with the cross 
sections obtained here for the bosonic lithium dimer 
( 7 Li 2 ). 



TABLE II: Elastic and total inelastic cross sections and rate 
coefficients for 7 Li + 7 Li2(vi,rn) at a collision energy of 0.928 
nK for different initial states of the molecule. 



D. Cold collisions 

At collision energies above about 1 mK, partial waves 
with end-over-end angular momentum I > contribute 
significantly to the overall cross sections. The conver- 
gence of the elastic and inelastic probabilities ak 2 /ir 
with total angular momentum J is shown in Figures 1151 
and [1751 for collisions with molecules in the (1,0) initial 
state. Partial waves J = to 10 suffice to converge 
cross sections up to ~ 580 mK, where the contribution of 
J = 10 to elastic and inelastic cross sections is 2.35% and 
0.0655%, respectively. It may be seen that, for each en- 
ergy, the inelastic cross sections are close to 1 for partial 
waves up to an energy-dependent cutoff value. The corre- 
sponding elastic cross sections oscillate around 1 until the 
cutoff is reached. Above 100 mK, convergence is slower 
for elastic probabilities, while below that it is slower for 
inelastic probabilities. 

The cross sections for partial waves with J > also 
show a simple energy dependence. In the Wigner regime 
the inelastic cross sections increase with slope I — 1/2 on 
a log-log scale, with / = J for rii = 0. At energies above 
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at an atom-diatom separation of 



FIG. 15: Elastic probabilities, ak 2 /-k, as a function of total 
angular momentum J for 7 Li + 7 Li2(v = l,n = 0). 




FIG. 16: Inelastic probabilities, ak 2 /it, as a function of total 
angular momentum J for 7 Li + 7 Li2(« = l,n = 0). 



the centrifugal barrier they decrease with slope —1. 

When several partial waves contribute significantly to 
the cross sections, the inelastic cross sections can be de- 
scribed semiquantitatively by classical Langevin capture 
theory. The assumption underlying this model is that 
the only condition for an inelastic event to take place 
is for the projectile to surmount the barrier of the effec- 
tive potential resulting from the repulsive centrifugal and 
attractive long-range potential. For a long-range interac- 
tion of the form —C s /R~ 8 the effective potential is 



KM = 



h 2 l{l + i) 



R s 



This has a barrier height of 



V = 

y max 



h 2 l(l + i) 



s/(s-2) 



(17) 



h 2 i{i 



i) 



l/( s -2) 



(19) 



Taking s = 6, C§ = 3085.54 £^Oq, and fi = 2mu/3, we 
obtain barrier heights of 2.78 mK, 14.4 mK, 40.8 mK, 
87.9 mK, 161 mK for I = 1 to 5, at atom-molecule sep- 
arations ranging from 94.25 ao for I = 1 to 47.89 ao for 
1 = 5. At these distances the interaction is dominated by 
dispersion and is nearly isotropic. The Langevin expres- 
sion for the inelastic cross section is 



O-inel(-E) 



S7T 

T 



(8-2)/ 8 



2/.s 



(20) 



The total inelastic cross sections are compared with this 
result in Figure [17] and it may be seen that there is in- 
creasingly good agreement between the full calculations 
and the Langevin model above about 30 mK. 

To test that the agreement between the Langevin 
model and our results is not fortuitous, we recalculated 
the cross sections using partial waves J = to 10 us- 
ing a different potential surface, designated ACVTZ. The 
ACVTZ potential was obtained in the same manner as 
described above, but with electronic energies calculated 
using a considerably smaller basis set (aug-cc-pCVTZ). 
The potential is significantly different, with a global D^h 
minimum 3873.37 cm -1 deep at r — 3.125 A and a 
minimum 930.29 cm -1 deep at r = 3.780 A. The two po- 
tentials have the same long-range form, joined using the 
switching function of Eq. ([7]). The elastic and inelastic 
cross sections calculated using the ACVTZ potential are 
included in Figure 1171 Once again the agreement with 
the Langevin model is excellent above ~ 30 mK, though 
the two potentials give cross sections that differ by about 
an order of magnitude in the Wigner regime. We found 
similar agreement with the model for collisions involving 
other initial molecular states. The Langevin cross sec- 
tions are independent of the initial state, particle masses 
and the details of short-range potential. 

Elastic cross sections are comparable to inelastic ones 
in the millikelvin region and can be as high as twice the 
inelastic ones. However, the ratio is never high enough to 
be favorable for the prospects of evaporative or sympa- 
thetic cooling. Whenever the Langevin model provides 
good estimates of inelastic cross sections, it implies that 
the inelastic probabilities are close to 1, which in turn 
implies that elastic cross sections (Eq. (jTU)) ) and inelastic 
cross sections (Eq. (fTTj) ') are not drastically different. 

The reaction 7 Li + 7 Li2 has no barrier for either lin- 
ear or perpendicular approaches. The deep well in the 
potential energy surface suggests that the reaction pro- 
ceeds via an insertion mechanism involving complex for- 
mation. It is therefore interesting to compare the scat- 
tering results obtained here with those for insertion re- 
actions that have been studied at ordinary temperatures 
(~ 100 meV) i 67 ' 68 i 69 ' 70 These reactions are also charac- 
terized by high inelastic and reactive probabilities. The 
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FIG. 17: [Color online] Elastic and total inelastic cross sec- 
tions for 7 Li + r hi2(vi = I, rii =0) on the uncontracted cc- 
pV5Z and aug-cc-pCVTZ basis set potentials and the inelastic 
cross sections in the Langevin model. 



final-state vibrational distributions decrease with v and 
the rotational distributions peak at high n for each final 
v. If all product states are equally probable, the product 
state distribution is proportional to the density of avail- 
able states, which depends on the rotational degeneracy 
and the density of translational states. The resulting 
differential cross sections display forward-backward sym- 
metry. This can be explained by formation of a collision 
complex whose decay is statistical^ It has been shown 
that the exact quantum scattering results for insertion 
reactions are in excellent agreement with results based 
on a quantum statistical theory at thermal energies^ 

The success of the Langevin model for total inelastic 
cross sections makes it worthwhile to investigate whether 
the product state distributions are well described by sta- 
tistical theory. The statistical prior distributions are 
given by 



Pv 



(2n + l)^E-E m , 



(21) 



where v and n refer to the final states and E is the total 
energy in the system. 

Figure [T8l compares the statistical predictions with full 
dynamics results in both the ultracold regime and at 580 
mK for 7 Li + 7 Li2 collisions for initial (v, n) — (3, 0). The 
ultracold results are quite different from the statistical 
predictions and show oscillations with n for each v simi- 
lar to those in Fig. [14] However, the oscillations start to 
wash out at higher energies and the rotational distribu- 
tions start to resemble the statistical ones, though they 
still favor low rotational quantum numbers. The vibra- 
tional distributions are clearly some way from statistical 
even at 580 mK and favor high v (smaller changes in 
vibrational quantum number). 

The corresponding results for differential cross sections 
are shown in Fig. 1191 Statistical theory predicts forward- 
backward symmetry, but this has not fully developed 
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FIG. 18: Final rotational distributions for 7 Li + 7 Li2(vi = 
3, rii — 0). Top panel: statistical prediction; center panel: 
ultracold regime (0.928 mK); bottom panel: collision energy 
of 580 mK. 



even at 580 mK. Similar results have been obtained in 
the study of cold collisions of K + K2^ The symmetry 
predicted by statistical theory arises due to phase can- 
cellations between cross terms in the sums over partial 
waves in Eq. ([8]) , and is incomplete at energies where only 
a few partial waves contribute. We expect that at higher 
collision energies, where more partial waves contribute, 
the symmetry observed in other deep-well systems will 
be recovered. 



E. Potential sensitivity 

Cross sections for elastic atom-atom collisions in the 
limit of zero collision energy depend on a single pa- 
rameter, the scattering length. The scattering length 
is very sensitive to details of the interaction potential 
and its properties and dependence on potential are well 
established.— In particular, the scattering length passes 
through a pole whenever there is a bound state of the 
diatomic potential curve at exactly zero energy, and the 
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FIG. 19: Center-of-mass differential cross sections for 7 Li + 
7 Li2(vi — 3, rii — 0) at a collision energy of 580 mK. Dif- 
ferential cross sections are integrated through the azimuthal 
angle and summed over the final states in each vibrational 
manifold. 



elastic cross section shows a corresponding peak of height 
4ir/k 2 (or 87r/fc 2 for identical bosons). 

A quantitative theoretical prediction of the scatter- 
ing length is within reach today only for the lightest 
diatomic systems and with the best available ab initio 
potentialsj 87 i 88 ' 89 For Li2, the sensitivity of scattering 
length to small differences between potentials is illus- 
trated in Table [T] For heavier systems the variation of 
scattering length with potential parameters is usually too 
fast to allow quantitative predictions until experimental 
data are available to help determine the scattering length. 

The absolute accuracy of the cross sections presented 
in this paper is limited by the quality of potential en- 
ergy surface. The accuracy of the electronic energies is 
degraded by the limitations of the correlation treatment 
and the size of the basis set used. Moreover, the effects 
of conical intersections and the influence of the excited 
electronic state have been neglected in our calculations. 
It is thus very worthwhile to investigate how changes in 
the potential energy surface affect the ultracold scatter- 
ing results. 

Since the potential for Li2 is accurately known, we 
modified the Li3 potential by scaling the nonadditive 
part of the potential by a multiplicative parameter A. 
The depth of the global minimum V m i n varies with the 
scaling factor as AV m i n /AX w 5150 cm -1 over the range 
studied. 

The dependence of the elastic cross sections on the 
scaling factor A for 7 Li + 7 Li2(vi = 0,nj = 0) at ~ 1 
nK is shown in Figure 1201 Whenever there is a bound 
state at exactly zero energy, the scattering length passes 
through a pole and the cross section shows a very high 
peak (cr w Air/k 2 ). Near each pole is a point where Su = 
1 so a = 0. It may be seen that the cross section passes 
through about 10 cycles for a 1% change in A. 




FIG. 20: Dependence of the elastic cross sections for 7 Li + 
7 Li2{vi — 0, rii — 0) on the scaling factor A of the nonadditive 
part of the potential. 



The mean value of the cross section can be estimated 
using the formula for the mean scattering length a given 
by Gribakin and Flambaum.— Taking the isotropic atom- 
molecule dispersion coefficient C§ = 3085.54 -E^a^, we 
obtain a = 5.83 x 10~ 13 cm 2 , which underestimates our 
accurate quantum results but gives the correct order of 
magnitude. 

For collisions of molecules in vibrationally excited 
states {vi,rii = 0) with Vi = 1 to 3, the dependence 
of elastic and inelastic cross sections at ~ 1 nK on A is 
shown in Figures [51] and [221 The dependence is again 
oscillatory, but the frequency and amplitude of the oscil- 
lations are smaller for higher initial vibrational quantum 
numbers Vi. The elastic (inelastic) cross sections vary 
across the range by 186% (136%) for v.^ = 1 and 64% 
(73%) for Vi — 3. The sharp peaks in cross sections asso- 
ciated with poles in scattering lengths are entirely washed 
out for higher values of i>; . 

The lack of strong resonant peaks as a function of 
well depth for Li + Li2 is an example of a general ef- 
fect recently discussed by Hutson. 90 Similar effects have 
recently been observed when tuning through zero-energy 
Feshbach resonances in He + NH using magnetic fields^ 
In the presence of inelastic scattering, the complex scat- 
tering length a = a — i/3 exhibits an oscillation rather 
than a pole when a bound-state crosses threshold. The 
amplitude of the oscillation is characterized by a reso- 
nant scattering length a rcs . If the bound state is coupled 
with comparable strengths to the elastic and inelastic 
channels, the oscillation may be of relatively low am- 
plitude and the peaks in cross sections are suppressed. 
This is the case in Li + Li 2 , where there are strong cou- 
plings between all the channels involved. For m = 1 there 
is enough resonant coupling to the elastic channel that 
peaks are still observed, but even these have an ampli- 
tude of only about a factor of 10. For higher Vi, there are 
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FIG. 21: Dependence of the elastic cross sections for 7 Li + 
7 hi2 (vi,rii = 0) on the scaling factor A of the nonadditive 
part of the potential. 
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FIG. 22: Dependence of the total inelastic cross sections for 
7 Li+ 7 Li2(«i, Hi = 0) on the scaling factor A of the nonadditive 
part of the potential. 



larger numbers of inelastic channels, each of which com- 
petes with the elastic channel and the resonant structures 
are progressively washed out. 

The suppression of cross-section peaks observed here 
for Li + hi2 contrasts with the situation commonly found 
for atom-atom scattering^ and for atom-molecule reac- 
tive scattering in systems such as F + For atom- 
atom collisions involving atoms in spin-stretched states, 
the resonant state is coupled to inelastic channels only 
by very week magnetic dipole-dipole interactions, but is 
coupled to the elastic channel by much stronger spin- 
exchange terms. This gives a large value of a Tes , corre- 



sponding to pole-like behavior of the scattering length 
and producing very large peaks in cross sections. Sim- 
ilarly, in F + H 2 the resonances studied^ are due to 
quasibound states in the entrance channel, while the "in- 
elastic" channels are actually reactive channels that are 
separated from the entrance well by a high barrier. This 
again gives very large peaks in cross sections. It is evi- 
dent that the factor that suppresses the resonant peaks 
in Li + Li2 is strong coupling to the inelastic channels. 

A sensitivity study similar to the one we report here 
was performed on the Na + Na2 system by Quemener 
et al^- They also observed oscillations that decreased in 
amplitude with initial v. 



F. Isotope effects 

A strong suppression of inelastic rates has been ob- 
served experimentally in fermionic systems tuned to a 
large and positive atom-atom scattering length^ 9 - ' 10 ! 93 
The molecular cloud consisting of weakly bound dimers 
of fermions exhibits remarkable stability against colli- 
sional decay. By contrast, the vibrational relaxation for 
weakly bound bosonic dimers is fast4&£ ' 18 i 94 The sup- 
pression of inelastic scattering for fermionic systems has 
been interpreted in terms of the requirements of Fermi 
statistics 

We have repeated all the above scattering calculations 
involving molecules in low-lying rovibrational states for 
the fermionic system 6 Li + 6 Li2. A brief account of the 
results was given in ref. [25|. For 6 Li + 6 Li2, the domi- 
nant contribution to cross sections in the ultracold regime 
comes from total angular momentum J = X~ , which is 
the partial wave that contains a contribution from U = 0. 
Cross sections at a collision energy of s=s 1 nK are shown 
in Table ITTT1 The individual calculated values are signif- 
icantly different from the bosonic case, but the effect is 
not due to suppression by Fermi statistics. Instead, the 
dominant effects is that the change in nuclear mass alters 
the scattering by a mechanism similar to a change in the 
potential surface, as discussed in the previous subsection. 
The cross sections calculated for individual initial states 
are thus once again scattered essentially randomly about 
a mean value. However, when we compare the results for 
a variety of initial rovibrational states, we see that there 
are no systematic differences between the overall rates 
for bosonic and fermionic systems for the low vibrational 
states that we have studied. 

Unusually low inelastic rates (about an order of magni- 
tude lower, with small Im(a)) were obtained for systems 
that have a small number of energetically accessible in- 
elastic channels, such as 6 Li + 6 Li2(t>i = 0,rti = 3) and 
7 Li + 7 Li 2 (fi = 0,rii = 2). 

The energy dependence of the cross sections in 
fermionic systems is qualitatively similar to that de- 
scribed above for bosonic systems. The Wigner regime 
is reached at energies around 10~ 5 K. For non-zero ini- 
tial end-over-end angular momenta, the simple relation- 
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TABLE III: Elastic and total inelastic cross sections and rate 
coefficients for 6 Li + 6 Li2(ui, rn) at a collision energy of 0.928 
nK for different initial states of the molecule. 
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FIG. 23: Final rotational distributions for 6 Li + 6 Li2(vi 
1, m = 1) at a collision energy in the Wigner regime. 
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ship of Eq. [16] holds at energies below the corresponding 
centrifugal barriers. At higher energies, the energy de- 
pendence of inelastic cross sections follows the Langevin 
model, Eq. [IB 

For spin-stretched 6 Li + 6 Li 2 collisions, only odd ro- 
tational levels of Li 2 are allowed. Figure [^3] shows the 
product rotational distribution for Vi = 1, = 1 in the 
Wigner regime, and it may be seen that the distribution 
is again oscillatory. At higher collision energies the oscil- 
lations are washed out in the sum over partial waves, as 
shown in Figure [24j However, as for the bosonic system, 
the rotational distribution is significantly non-statistical. 
The forward-backward symmetry of the differential cross 
sections predicted by statistical theory is also not present 
for 6 Li + 6 Li2 at 116 mK, as seen in Figure [25] 



G. Reactions in isotopic mixtures 

A greater variety of atom-diatom collisions may be 
observed in mixtures. Isotopic mixtures of 6 Li and 7 Li 
are of considerable interest , 95 i 96 i 97 and one can envisage 
creating either homonuclear or heteronuclear dimers in 



FIG. 24: Final rotational distributions for 6 Li + 6 Li2{vi 
l,n, = 1) at a collision energy of 116 mK. 



such a mixture. The atom-diatom collision systems of 
interest are 



7 Li- 



'Li 



6 Li 7 Li, 



5 Li 2 . 



6 Li + 7 Li 2 , 



6 Li + 6 Li 7 Li, 



(22) 
(23) 
(24) 
(25) 



We have studied all four of the systems |22|) to (|25p and 
gave a preliminary report in ref. 26. The novel feature 
in these systems is the possibility of a chemical reaction 
in which the reactants and products are distinguishable. 
The low- lying energy levels of 7 Li 2 , 6 Li 7 Li and 6 Li 2 are 
shown relative to the Li 2 potential minimum in Fig. 1261 
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FIG. 25: Center-of-mass differential cross sections for 6 Li + 
6 Li2(vi = 1, m = 1) at a collision energy of 116 mK. Differen- 
tial cross sections are integrated through the azimuthal angle 
and summed over the final states in each vibrational manifold 
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FIG. 26: The low-lying energy rotational levels of 7 Li2, 6 Li 7 Li 
and 6 Li2 for v = 0, relative to the Li2 potential minimum. 
Only levels with even n are shown for 7 Li2 and only levels 
with odd n are shown for 6 Li2. 



Because of the differences in zero-point energy, the reac- 
tions are exoergic even for ground-state molecules for the 
two systems involving atomic 7 Li, (|22p and (|23p. 

For 6 Li 7 Li, both even and odd rotational levels are al- 
lowed even for spin-stretched states. The higher density 
of states makes these systems harder to treat computa- 
tionally. Our calculations on these systems were therefore 
restricted to the J = and 1 partial waves and are thus 
converged only at ultracold collision energies. 

We will consider the systems that are reactive for 
ground-state molecules first. The cross sections, rate 
coefficients and scattering lengths for collisions of 7 Li 
with 6 Li 7 Li in several different initial molecular states 



are summarized in Table IIVI For molecules in excited 
vibrational or rotational states, the cross sections are 
comparable to those found for the homonuclear systems. 
However, for 6 Li 7 Li(w = Q,n = 0) the cross section for 
reaction is a factor of about 50 smaller than for other 
states. This corresponds to a reactive rate coefficient 
of 4.05 x 10~ 12 cm 3 s^ 1 and Im(a) = 0.233 A. In this 
case there is only a single exoergic channel, with a ki- 
netic energy release of 1.822 K. The small cross section 
is probably due to the small volume of available phase 
space for reaction. 

The increased density of rotational states for 6 Li 7 Li 
produces an increased density of Feshbach resonances. 
The energy-dependence of the eigenphase sum for 7 Li + 
6 Li 7 Li is shown in the top panel of Figure [571 There is a 
Feshbach resonance clearly visible at 225 mK and another 
at about 470 mK, though the latter overlaps the thresh- 
old for opening the n — 2 reactive channel at 477.2 mK. 
There are further overlapping resonances near 630 mK 
and 725 mK. The corresponding structures in the elastic 
and inelastic s-wave cross sections are shown in the lower 
two panels of Figure [271 The elastic cross section shows 
a minimum near 120 mK, where the phase shift obtained 
from the elastic S-matrix element passes through a mul- 
tiple of 7r. It dips to a very small value near the Feshbach 
resonance at 225 mK, but not near the resonance near 470 
mK. The reactive cross section, by contrast, dips close to 
zero at both Feshbach resonances. This suggests the in- 
teresting possibility of reducing the inelastic/elastic cross 
section ratio by tuning close to a Feshbach resonance. 

Reactive collisions of 7 Li with 6 Li2 are more strongly 
exoergic because for spin-stretched states the 6 Li2 
molecule is restricted to odd rotational levels while the 
6 Li 7 Li molecule is not. Because of this collisions of 7 Li 
with 6 Li2 in its n = 1 ground state can produce 6 Li 7 Li 
in n = 0, 1 or 2. The calculated cross sections, rate coef- 
ficients and scattering lengths are summarized in Table 
fVl The reactive cross sections for 6 Li2 (n = 1) are again 
reduced compared to those for other states, but not as 
dramatically as for 7 Li + 6 Li 7 Li (n = 0). 

Ultracold collisions of 6 Li 7 Li and 7 Li2 with atomic 6 Li 
are nonreactive unless the molecule is initially in a rovi- 
brationally excited state. This is because for these sys- 
tems the zero-point energy is greater in the products than 
in the reactants. The calculated cross sections, rate coef- 
ficients and scattering lengths for these systems are sum- 
marized in Tables IVTl and IVTIl The patterns of cross sec- 
tions for rovibrationally excited initial states are fairly 
similar to those for the homonuclear systems, with no 
particularly small inelastic or reactive rates. The novel 
feature of these systems really lies in the fact that the 
reactions are forbidden for ground-state molecules, so 
that sympathetic cooling of molecules by contact with 6 Li 
atoms might be envisaged. For example, if 6 Li2 molecules 
are produced in the presence of both 6 Li and 7 Li, they 
can react to form (relatively hot) 6 Li 7 Li. If the 7 Li atoms 
are then removed, the molecules might be cooled by elas- 
tic collisions with 6 Li. 
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Vi, n t cr c i as (cm 2 ) q-jnei (cm 2 ) g- roac (cm 2 ) fci oss (cm 3 s -1 ) Re(a) (nm) -Im(a) (11m) 

0, 1.77 x 1(T 13 - 2.20 x lO -11 4.05 x 10~ 12 1.19 0.0233 

1, 1.32 x 10~ 12 8.79 x 10~ 10 2.56 x 10~ 1() 2.09 x 10~ 10 3.01 1.19 

2, 1.32 x 10~ 12 1.37 x 10~ 9 1.02 x 10~ 9 4.40 x 10~ 10 2.04 2.52 

3, 1.09 x 10~ 12 1.30 x 10~ 9 8.74 x 10~ lu 4.00 x 10~ 10 L86 2.29 

0, 1 4.81 x 10" 12 4.30 x 10" 10 2.32 x 10" 9 5.06 x 10~ 10 5.46 2.90 

1, 1 1.52 x 10~ 12 1.85 x 10~ 9 9.47 x 10~ 10 5.14 x 10~ 10 1.85 2.94 

2, 1 9.09 x 10~ 13 1.22 x 10~ 9 4.45 x 10~ 10 2.93 x 10~ 10 2.10 1.68 

3, 1 9.51 x 10~ 13 1.67 x 10~ 9 6.28 x 10~ 10 4.23 x 10~ 10 1.71 2.15 



TABLE IV: Cross sections and related parameters for 7 Li + 6 Li 7 Li (vi,ni) at a collision energy of 0.928 nK for different initial 
states of the molecule. 



Vj, rij cr c i as (cm 2 ) o-jnei (cm 2 ) q- rca c (cm 2 ) fci oss (cm 3 s -1 ) Re(q) (nm) -Im(a) (nm) 

0, 1 1.17 x 10~ 12 2.34 x 10~ 10 4.37 x 10~ n 3.04 0.244 

1, 1 1.62 x 10~ 12 1.00 x 10~ 9 1.76 x 10~ 9 5.15 x 10~ 10 2.17 2.86 

2, 1 6.61 x 10~ 13 2.66 x 10~ 10 1.12 x 10~ 9 2.58 x 10~ 10 1.41 1.81 

3, 1 8.05 x 10~ 13 5.29 x 10~ 10 1.07 x 10~ 9 2.98 x 10~ 10 1.91 1.66 



TABLE V: Cross sections and related parameters for 7 Li + 6 Li2 (vi,rn) at a collision energy of 0.928 nK for different initial 
states of the molecule. 



Vi, rij cr c i as (cm 2 ) o-jnei (cm 2 ) Uroac (cm 2 ) fci oss (cm 3 s -1 ) Re(q) (nm) -Im(a) (nm) 

0, 4.71 x 10~ 12 6.12 

1, 6.43 x 10~ 13 1.15 x 10~ 9 2.12 x 10 -10 2.64 x 10~ 10 1.81 1.36 

2, 1.12 x 10~ 12 1.39 x 10~ 9 4.37 x 10~ 10 3.53 x 10~ 10 2.36 1.82 

3, 1.54 x 10~ 12 1.34 x 10~ 9 9.49 x 10~ lu 4.43 x 10~ 10 2M 2.29 

0, 1 2.74 x 10" 12 5.84 x 10" 10 1.13 x 10~ 10 4.64 0.586 

1, 1 7.56 x 10~ 13 2.04 x 10~ 9 2.72 x 10~ 10 4.48 x 10~ 10 0.828 2.31 

2, 1 1.19 x 10~ 12 2.07 x 10~ 9 7.21 x 10~ 10 5.40 x 10~ 10 1.32 2.79 

3, 1 1.27 x 10~ 12 1.55 x 10~ 9 7.76 x 10~ 10 4.51 x 10~ 10 2.16 2.33 



TABLE VI: Cross sections and related parameters for 6 Li + 6 Li 7 Li (vi,rii) at a collision energy of 0.928 nK for different initial 
states of the molecule. 



Vi, n t <7 c i as (cm 2 ) o-jnei (cm 2 ) q- r eac (cm 2 ) fci oss (cm 3 s -1 ) Re(q) (nm) -Im(a) (nm) 

0, 1.29 x 10~ 12 3.20 

1, 1.33 x 10~ 12 5.70 x 10~ 10 8.85 x 10~ 10 2.79 x 10~ 10 2.90 1.47 

2, 1.14 x 10~ 12 1.39 x 10~ 9 1.37 x 10~ 9 5.28 x 10~ 10 1.16 2.79 

3, 1.51 x 10~ 12 1.00 x 10~ 9 1.42 x 10~ 9 4.63 x 10~ 10 2.45 2.45 



TABLE VII: Cross sections and related parameters for 6 Li + 7 Li2 (vi,ni) at a collision energy of 0.928 nK for different initial 
states of the molecule. 
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FIG. 27: Energy dependence of eigenphase sum (top panel), 
elastic cross section (center panel) and inelastic/reactive cross 
section (bottom panel) for collisions of 6 Li 7 Li (v = 0, n = 0) 
with 7 Li. 



IV. CONCLUSIONS 



Quantum dynamics calculations were carried out in a 
time-independent reactive scattering formalism based on 
an expansion in pscudohypcrsphcrical harmonics. This 
approach makes it straightforward to include the required 
symmetry with respect to exchange of identical bosons 
( 7 Li atoms) or fermions ( 6 Li atoms). For the homonu- 
clear systems 7 Li + 7 Li 2 and 6 Li + 6 Li 2 , the products of 
reactive and non-reactive collisions are indistinguishable. 
Nevertheless, all collisions involving vibrationally excited 
Li 2 result in very fast vibrational relaxation (quench- 
ing). For the low- lying vibrationally excited states stud- 
ied here, there is no systematic suppression of quenching 
rates for collisions involving fermion dimers ( 6 Li 2 ). This 
contrasts with the situation found both experimentally 
and theoretically for fermion dimers produced by Fesh- 
bach resonance tuning in their highest vibrational state. 

For collisions involving mixtures of Li isotopes, reactive 
and nonreactive collision outcomes can be distinguished. 
In particular, for collisions of 7 Li atoms with either 6 Li 2 
or 6 Li 7 Li, exoergic reactive collisions are possible even 
for vanishing kinetic energy because of the difference in 
zero-point energy between reactants and products. These 
reactive processes are generally very fast, except for the 
case of 7 Li + 6 Li 7 Li (v = 0, n = 0) at collision ener- 
gies below 477 mK, for which there is only one reactive 
channel. 

The Li + Li 2 collision systems have a rich structure of 
scattering resonances at low energies. These are princi- 
pally rotational Feshbach resonances due to quasibound 
states that correlate with rotationally excited diatomic 
molecules. Tuning such Feshbach resonances with ap- 
plied magnetic fields offers possibilities for controlling ul- 
tracold molecular collisions, and this will be a topic for 
future work. 



We have developed a new potential energy surface for 
quartet Li + Li 2 from high-level electronic structure cal- 
culations and used it to carry out quantum dynamics cal- 
culations on elastic, inelastic and reactive collisions on Li 
+ Li 2 under cold and ultracold collisions. The potential 
energy surface was calculated using RCCSD(T) calcula- 
tions with all electrons correlated. The surface was cal- 
culated on a grid in pure bond- length coordinates. It was 
interpolated at short range using the interpolant-moving 
least-squares method and merged with a long-range form 
that is correct at both the atom-diatom and the atom- 
atom-atom dissociation limits. 

The potential energy surface has a deep triangular min- 
imum that lies far below the atom-diatom energy and al- 
lows barrierless atom exchange reactions. The well depth 
is a factor of 4 deeper than would be expected on the basis 
of a pairwise-additive sum of triplet Li 2 potentials. The 
surface exhibits a seam of conical intersections at linear 
geometries, at energies close to that at three separated 
Li atoms. 



We have investigated the sensitivity of clastic and in- 
clastic cross sections to variations in the interaction po- 
tential. Such variations tune Feshbach resonances across 
thresholds, and our expectation was that the cross sec- 
tions would show very large peaks when this occurred. 
This is indeed the case for collisions of ground-state 
molecules where only elastic scattering is possible. For 
such collisions the elastic cross sections increase by many 
orders of magnitude close to each resonance, and the den- 
sity of states is such that a 1% change in the potential 
depth sweeps about 10 resonances across threshold. How- 
ever, for collisions where inelastic scattering is also pos- 
sible the resonant behavior is strongly suppressed. The 
amplitude of the oscillations in cross sections decreases 
dramatically with increasing initial vibrational quantum 
number, and for v = 3 the cross section changes by no 
more than a factor of 2 as resonances cross threshold. 
This shows that, in some cases, cross sections are much 
less sensitive to small variations in the potential than was 
previously expected. 



20 



Acknowledgments 

MTC is grateful for sponsorship from the University of 
Durham and Universities UK. 



1 J. M. Hutson and P. Soldan, Int. Rev. Phys. Chem. 25, 
497 (2006). 

2 K. M. Jones, E. Tiesinga, P. D. Lett, and P. S. Julienne, 
Rev. Mod. Phys. 78, 483 (2006). 

3 T. Kohler, K. Goral, and P. S. Julienne, Rev. Mod. Phys. 
78, 1311 (2006). 

4 E. A. Donley, N. R. Claussen, S. T. Thompson, and C. E. 
Wieman, Nature 417, 529 (2002). 

5 J. Herbig, T. Kraemer, M. Mark, T. Weber, C. Chin, H. C. 
Nagerl, and R. Grimm, Science 301, 1510 (2003). 

6 K. Xu, T. Mukaiyama, J. R. Abo-Shaeer, J. K. Chin, 
D. E. Miller, and W. Ketterle, Phys. Rev. Lett. 91, 210402 
(2003). 

7 S. Diirr, T. Volz, A. Marte, and G. Rempe, Phys. Rev. 
Lett. 92, 020406 (2004). 

8 K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. 
Rev. Lett. 91, 080406 (2003). 

9 J. Cubizolles, T. Bourdel, S. J. J. M. F. Kokkelmans, G. V. 
Shlyapnikov, and C. Salomon, Phys. Rev. Lett. 91, 240401 
(2003). 

10 S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, 
C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. Lett. 
91, 240402 (2003). 

11 S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, 
S. Riedl, C. Chin, J. H. Denschlag, and R. Grimm, Sci- 
ence 302, 2101 (2003). 

12 M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 
Raupach, S. Gupta, Z. Hadzibabic, and W. Ketterle, Phys. 
Rev. Lett. 91, 250401 (2003). 

13 M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 
(2003). 

14 P. Staanum, S. D. Kraft, J. Lange, R. Wester, and M. Wei- 
demiiller, Phys. Rev. Lett. 96, 023201 (2006). 

15 J. M. Sage, S. Sainis, T. Bergeman, and D. DeMille, Phys. 
Rev. Lett. 94, 203001 (2005). 

16 J. M. Hutson and P. Soldan, Int. Rev. Phys. Chem. 26, 1 
(2007). 

17 R. Wynar, R. S. Freeland, D. J. Han, C. Ryu, and D. J. 
Heinzen, Science 287, 1016 (2000). 

18 T. Mukaiyama, J. R. Abo-Shaeer, K. Xu, J. K. Chin, and 
W. Ketterle, Phys. Rev. Lett. 92, 180402 (2004). 

19 N. Zahzam, T. Vogt, M. Mudrich, D. Comparat, and P. Fil- 
let, Phys. Rev. Lett. 96, 023202 (2006). 

20 P. Soldan, M. T. Cvitas, J. M. Hutson, P. Honvault, and 
J. M. Launay, Phys. Rev. Lett. 89, 153201 (2002). 

21 G. Quemener, P. Honvault, and J. M. Launay, Eur. Phys. 
J. D 30, 201 (2004). 

22 D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. 
Rev. Lett. 93, 090404 (2004). 

23 D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. 
Rev. A 71, 012708 (2005). 

24 G. Quemener, P. Honvault, J. M. Launay, P. Soldan, D. E. 
Potter, and J. M. Hutson, Phys. Rev. A 71, 032722 (2005). 

25 M. T. Cvitas, P. Soldan, J. M. Hutson, P. Honvault, and 
J. M. Launay, Phys. Rev. Lett. 94, 033201 (2005). 



26 M. T. Cvitas, P. Soldan, J. M. Hutson, P. Honvault, and 
J. M. Launay, Phys. Rev. Lett. 94, 200402 (2005). 

27 J. Higgins, T. Hollebeek, J. Reho, T. S. Ho, K. K. 
Lehmann, H. Rabitz, G. Scoles, and M. Gutowski, J. 
Chem. Phys. 112, 5751 (2000). 

28 P. Soldan, M. T. Cvitas, and J. M. Hutson, Phys. Rev. A 
67, 054702 (2003). 

29 D. A. Brue, X. Li, and G. A. Parker, J. Chem. Phys. 123, 
091101 (2005). 

30 H.-J. Werner and W. Meyer, J. Chem. Phys. 73, 2342 

(1980) . 

31 H.-J. Werner and W. Meyer, J. Chem. Phys. 74, 5794 

(1981) . 

32 H.-J. Werner and P. J. Knowles, J. Chem. Phys. 82, 5053 
(1985). 

33 P. J. Knowles and H.-J. Werner, Chem. Phys. Lett. 115, 
259 (1985). 

34 H.-J. Werner, Adv. Chem. Phys. 69, 1 (1987). 

35 H.-J. Werner, P. J. Knowles, R. Lindh, M. Schiitz, 
P. Celani, T. Korona, F. R. Manby, G. Rauhut, 
R. D. Amos, A. Bernhardsson, et al., Molpro, ver- 
sion 2002.6, a package of ab initio programs (2003), see 
http://www.molpro.net. 

36 Emsl gaussian basis set order form at 
http://www.emsl.pnl.gov/forms/basisform.html. 

37 H.-J. Werner and E. A. Reinsch, J. Chem. Phys. 76, 3144 

(1982) . 

38 P. J. Knowles and H.-J. Werner, Chem. Phys. Lett. 145, 
514 (1988). 

39 H.-J. Werner and P. J. Knowles, J. Chem. Phys. 89, 5803 
(1988). 

40 F. D. Colavecchia, J. P. Burke, W. J. Stevens, M. R. 
Salazar, G. A. Parker, and R. T Pack, J. Chem. Phys. 
118, 5484 (2003). 

41 T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989). 

42 M. A. Iron, M. Oren, and J. M. L. Martin, Mol. Phys. 101, 
1345 (2003). 

43 Z. C. Yan, J. F. Babb, A. Dalgarno, and G. W. F. Drake, 
Phys. Rev. A 54, 2824 (1996). 

44 R. W. Molof, H. L. Schwartz, T. M. Miller, and B. Beder- 
son, Phys. Rev. A 10, 1131 (1974). 

45 C. E. Moore, Atomic Energy Levels, vol. 1 (National Bu- 
reau of Standards, Washington, 1971). 

46 S. F. Boys and F. Bernardi, Mol. Phys. 19, 553 (1970). 

47 T. S. Ho and H. Rabitz, J. Chem. Phys. 104, 2584 (1996). 

48 T. S. Ho and H. Rabitz, J. Chem. Phys. 113, 3960 (2000). 

49 Z. C. Yan, A. Dalgarno, and J. F. Babb, Phys. Rev. A 55, 
2882 (1997). 

50 C. Linton, F. Martin, A. J. Ross, I. Russier, P. Crozet, 
A. Yiannopoulou, L. Li, and A. M. Lyyra, J. Mol. Spec- 
trosc. 196, 20 (1999). 

51 E. R. I. Abraham, W. I. McAlexander, C. A. Sackett, and 
R. G. Hulet, Phys. Rev. Lett. 74, 1315 (1995). 

52 T. J. Lee and P. R. Taylor, Int. J. Quantum Chem. pp. 
199-207 (1989). 



21 



53 M. T. Cvitas, P. Soldan, and J. M. Hutson, Mol. Phys. 
104, 23 (2006). 

54 T. Ishida and G. C. Schatz, J. Comput. Chem. 24, 1077 
(2003). 

55 T. Ishida and G. C. Schatz, J. Chem. Phys. 107, 3558 

(1997) . 

56 T. Ishida and G. C. Schatz, Chem. Phys. Lett. 298, 285 

(1998) . 

57 J. Ischtwan and M. A. Collins, J. Chem. Phys. 100, 8080 
(1994). 

58 B. M. Axilrod and E. Teller, J. Chem. Phys. 11, 299 
(1943). 

59 Y. Muto, Proc. Phys. Math. Soc. Japan 17, 629 (1943). 

60 R. J. Bell, J. Phys. B - At. Mol. Opt. Phys. 3, 751 (1970). 

61 M. B. Doran and I. J. Zucker, J. Phys. C - Solid St. Phys. 
4, 307 (1971). 

62 W. L. Bade, J. Chem. Phys. 27, 1280 (1957). 

63 W. L. Bade, J. Chem. Phys. 28, 282 (1958). 

64 K. T. Tang and J. P. Toennies, J. Chem. Phys. 80, 3726 
(1984). 

65 S. H. Patil and K. T. Tang, J. Chem. Phys. 106, 2298 
(1997). 

66 M. Rerat and B. Bussery-Honvault, Mol. Phys. 101, 373 
(2003). 

67 P. Honvault and J. M. Launay, J. Chem. Phys. Ill, 6665 

(1999) . 

68 P. Honvault and J. M. Launay, J. Chem. Phys. 114, 1057 
(2001). 

69 F. J. Aoiz, L. Banares, J. F. Castillo, M. Brouard, W. Den- 
zer, C. Vallance, P. Honvault, J. M. Launay, A. J. Dobbyn, 
and P. J. Knowles, Phys. Rev. Lett. 86, 1729 (2001). 

70 P. Honvault and J. M. Launay, Chem. Phys. Lett. 370, 
371 (2003). 

71 J. M. Launay and M. LeDourneuf, Chem. Phys. Lett. 163, 
178 (1989). 

72 D. E. Manolopoulos, J. Chem. Phys. 85, 6425 (1986). 

73 N. Balakrishnan, V. Kharchenko, R. C. Forrey, and A. Dal- 
garno, Chem. Phys. Lett. 280, 5 (1997). 

74 E. P. Wigner, 73, 1002 (1948). 

75 C. J. Ashton, M. S. Child, and J. M. Hutson, J. Chem. 
Phys. 78, 4025 (1983). 

76 J. M. Hutson, C. J. Ashton, and R. J. Le Roy, J. Phys. 
Chem. 87, 2713 (1983). 

77 B. Stewart, P. D. Magill, T. P. Scott, J. Derouard, and 



D. E. Pritchard, Phys. Rev. Lett. 60, 282 (1988). 

78 P. D. Magill, B. Stewart, N. Smith, and D. E. Pritchard, 
Phys. Rev. Lett. 60, 1943 (1988). 

79 W. J. Hoving and R. Parson, Chem. Phys. Lett. 158, 222 
(1989). 

80 R. C. Forrey, Phys. Rev. A 66, 023411 (2002). 

81 P. M. Florian, M. Hoster, and R. C. Forrey, Phys. Rev. A 
70, 032709 (2004). 

82 K. Tilford, M. Hoster, P. M. Florian, and R. C. Forrey, 
Phys. Rev. A 69, 052705 (2004). 

83 R. C. Forrey, N. Balakrishnan, A. Dalgarno, M. R. Hag- 
gerty, and E. J. Heller, Phys. Rev. Lett 82, 2657 (1999). 

84 R. D. Levine and R. B. Bernstein, Molecular Reaction Dy- 
namics and Chemical Reactivity (Oxford University Press, 
1987). 

85 E. J. Rackham, T. Gonzalez-Lezana, and D. E. 
Manolopoulos, J. Chem. Phys. 119, 12895 (2003). 

86 G. F. Gribakin and V. V. Flambaum, Phys. Rev. A 48, 
546 (1993). 

87 F. X. Gadea, T. Leininger, and A. S. Dickinson, Eur. Phys. 
J. D 15, 251 (2001). 

88 F. X. Gadea, T. Leininger, and A. S. Dickinson, J. Chem. 
Phys. 117, 7122 (2002). 

89 A. S. Dickinson, F. X. Gadea, and T. Leininger, J. Phys. 
B - At. Mol. Opt. Phys. 37, 587 (2004). 

90 J. M. Hutson, New J. Phys. 9, 152 (2007). 

91 M. L. Gonzalez-Martinez and J. M. Hutson, Phys. Rev. A 
75, 022702 (2007). 

92 E. Bodo, F. A. Gianturco, N. Balakrishnan, and A. Dal- 
garno, J. Phys. B - At. Mol. Opt. Phys. 37, 3641 (2004). 

93 C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 
92, 083201 (2004). 

94 J. Stenger, S. Inouye, M. R. Andrews, H. J. Miesner, D. M. 
Stamper-Kurn, and W. Ketterle, Phys. Rev. Lett. 82, 2422 
(1999). 

95 A. G. Truscott, K. E. Strecker, W. I. McAlexander, G. B. 
Partridge, and R. G. Hulet, Science 291, 2570 (2001). 

96 F. Schreck, G. Ferrari, K. L. Corwin, J. Cubizolles, 
L. Khaykovich, M. O. Mewes, and C. Salomon, Phys. Rev. 
A 6401, 011402 (2001). 

97 L. Khaykovich, F. Schreck, J. Cubizolles, T. Bourdel, K. L. 
Corwin, G. Ferrari, and C. Salomon, Physica B 329, 13 
(2003). 



